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Abstract 

Non-standard neutrino interactions (NSI) involved in neutrino propagation inside Earth matter 
could potentially alter atmospheric neutrino fluxes. In this work, we look at the impact of these 
NSI on the signal at the ICAL detector to be built at the India-based Neutrino Observatory (INO). 
We show how the sensitivity to the neutrino mass hierarchy of ICAL changes in the presence of 
NSI. The mass hierarchy sensitivity is shown to be rather sensitive to the NSI parameters e e/1 and 
e er , while the dependence on e /ir and e TT is seen to be very mild, once the y 2 is marginalised 
over oscillation and NSI parameters. If the NSI are large enough, the event spectrum at ICAL 
is expected to be altered and this can be used to discover new physics. We calculate the lower 
limit on NSI parameters above which ICAL could discover NSI at a given C.L. from 10 years of 
data. If NSI were too small, the null signal at ICAL can constrain the NSI parameters. We give 
upper limits on the NSI parameters at any given C.L. that one is expected to put from 10 years 
of running of ICAL. Finally, we give C.L. contours in the NSI parameter space that is expected to 
be still allowed from 10 years of running of the experiment. 
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I. INTRODUCTION 


The 50 kton magnetised iron detector (ICAL) to be built at the India-based neutrino 
observatory (INO) will be mainly observing muon neutrinos coming from the Earth’s atmo¬ 
sphere jl| . Amongst the most important goals of this experiment is the determination of 
the neutrino mass hierarchy (MH) through the observation of Earth matter effects in the 
expected data sample. For the mass-squared difference > 0 (normal hierarchy) 1 , one 

expects matter enhanced oscillations in the neutrino channel in the energy range around 
(5-10) GeV, while the antineutrino channel does not experience any such matter induced 
enhancement. On the other hand for A m\ v < 0 (inverted hierarchy), matter enhanced 
oscillations are expected in the antineutrino channel, while the neutrino channel does not 
obtain any such enhancement. ICAL being magnetised will be able to tell its neutrino signal 
from its antineutrino signal, giving the detector an added handle on measuring these Earth 
matter effects, and hence, the neutrino mass hierarchy. The sensitivity reach of this experi¬ 
ment for measuring standard neutrino oscillation parameters have been studied extensively 
in Refs. 

So far, there has been no signal of physics beyond the Standard Model in any of the 
accelerator-based experiments including LHC. However, we have unambiguous evidence from 
complimentary experiments that the Glashow-Weinberg-Salam model of elementary parti¬ 
cles is at best a low-energy effective theory and that there exists physics beyond the Standard 
Model. The Standard Model of particle physics is unable to convincingly explain data from 
neutrino oscillation experiments. In addition, it also fails to provide explanation of the exis¬ 
tence of dark matter and dark energy in the Universe, as well as baryogenesis. There are also 
theoretical issues with the Standard Model which demand its extension. Any extension of 
the model could include addition to its gauge or particle sector, or both. It is therefore per¬ 
tinent to envisage that such an extended theory would also have new (effective) interactions 
between the particles, beyond what is included in the Standard Model. Such interactions are 
expected to change the predicted outcome of experiments and existing data can be used to 
put limits on the strength of these interactions. In this work, we are primarily interested in 
non-standard interactions (NSI) affecting neutrino oscillations as neutrinos propagate inside 

1 We define Am?- = mf — m?. 
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Earth matter. These NSI, if present, would modify the transport of atmospheric neutrinos 
inside Earth matter, and hence alter the signal at the ICAL detector. 

The currently running Super-Kamiokande (SK) atmospheric neutrino experiment has 
looked for possible presence of these NSI in its event sample, and has found the data to 
be completely consistent with the Standard Model. Through a statistical analysis, the SK 


collaboration converts this into an upper bound on the relevant NSI parameters 


0 . 


Ex¬ 


pected constraints from other (future) atmospheric neutrino experiments have been studied 


previously in the literature, see e.g. Refs. js4l3j| (see Ref. 14-[l9j for earlier works). 


In this paper, we will study in detail the impact of NSI on the atmospheric neutrino signal 
in the ICAL detector at INO. We analyse the prospective data at ICAL in terms of the 
measured muon energy and muon angle through a binned y 2 analysis. ICAL is expected to 
also measure the energy deposited in the associated hadron shower. Inclusion of the hadron 
energy information improves the energy reconstruction of the events and hence in general 


improves the sensitivity of ICAL jf|. We have not included the hadron energy information 
in this work. This is being studied in a follow-up work by the INO collaboration 20]. We 
use the Nuance event generator with the ICAL detector geometry for generating muons 
from atmospheric neutrinos. The ICAL energy and angle resolutions and reconstruction 
and charge identification efficiencies are obtained from the Geant-based detector simulation 
code developed for ICAL. We generate muon events in the range (1-100) GeV and show the 
increase in the sensitivity to NSI parameters with the increase of the neutrino energy, and 
hence the muon energy, as was pointed out in Ref. 10]. We will quantify the extent of this 
modification in the expected muon signal at ICAL. Using that we will study the expected 
limits that ICAL could impose on NSI parameters if there is no evidence of NSI in the data. 
If on the other hand the NSI parameters are large enough, we would see a signal of new 
physics in the ICAL data. We give the lower limit on the NSI parameters which is needed 
in order to allow their discovery in ICAL at any given C.L. Likewise, the presence of NSI 
could change the sensitivity of ICAL to other neutrino oscillation parameters. In particular, 
we will show how the NSI parameters alter the mass hierarchy sensitivity in this class of 
experiments and present the revised sensitivity limits. 

The paper is organised as follows. In Section we discuss the neutrino oscillation 
probabilities in the presence of NSI. In particular, we study the impact of NSI on the 
difference in the probabilities between the NH and IH cases. In Section IIII1 we describe the 
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ICAL detector, our simulation techniques, and the statistical analysis procedure. We present 
our main results in Section IIVI All results are shown for 500 kton-year of data in ICAL. 
Subsection 11V Al is devoted to study the impact of NSI on the mass hierarchy sensitivity 
of ICAL. In Subsections IIV B1 II V Cl and IIV D1 we discuss the expected constraints on NSI 
parameters, the expected lower limit allowing for discovery of NSI, and the allowed areas in 
NSI parameters space, respectively. We end in Section lYl with our conclusions. 


II. IMPACT OF NSI ON OSCILLATION PROBABILITIES 


As outlined in the Introduction, an extension of the Standard Model of particle physics 
in its gauge sector and/or particle sector is likely to give rise to additional (effective) interac¬ 
tions between Standard Model particles. In particular, in this work, we are concerned with 


such interactions experienced by the neutrinos when_t 


This effective term in the Lagrangian is of the form 21 


hey 


24] 


propagate inside Earth matter. 


£ NS i = -2V2G F e f a c J^P L vp)(f^P c f ), 


( 1 ) 


where / is a fermion, Pc = (1 ±7 5 )/2 {C = R, L) are the chiral projection operators, Gp is 
the Fermi constant, and are the corresponding NSI parameters. Since Earth matter is 
made up of the first generation fermions only, the NSI parameters corresponding to e, u, and 
d are the only ones which contribute towards modifying the neutrino propagation inside the 
Earth. For the neutral-current NSI what is relevant is the sum e afj = e^ + e^. Furthermore, 
since only the incoherent sum of the NSI contributions is important, we combine the NSI 
effects coming from e^, e“^, and as 




£ 


"u 

n P af} 


( 2 ) 


where rif is the number density of the fermion / and we have normalised the effective contri¬ 
bution to the number density of electrons in Earth. 2 While, in principle, the NSI parameters 
are complex, we consider only values which have phases either 0 and 7r. Throughout this 


2 Note that the convention followed in defining the NSI parameters is crucial to interpret the actual con¬ 
straints on them from a given experiment and there are some places in the literature where a different 
convention has been followed (see Ref. |25| for a discussion). 
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work, we will use this assumption. Note that we sometimes refer to the NSI parameters as 
e a p for simplicity, even though we work with only the real values of these parameters. 

Each of the NSI parameters has been constrained from existing data. The corresponding 


model-independent upper bounds at 90 % C.L. are 

( 

l e adl < 


26] 


4.2 0.33 3.0 
0.33 0.068 0.33 
3.0 0.33 21 


\ 


( 3 ) 


where we have arranged the parameters in the form of a matrix with the rows and columns 
corresponding to {e, /i, r}. Note that only the NSI parameter e MA[ is well constrained in this 
phenomenological approach, while constraints on all other NSI parameters are rather loose. 
In particular, large values of e ee , e er , and e TT are still allowed. These bounds are generally 
referred to in the literature as indirect bounds as these bounds on parameters affecting 
neutrino oscillations come from non-neutrino experiments. The only neutrino experiment 
that has provided bounds on (some) of these parameters that are better than these indirect 
bounds is the SK experiment, which puts direct bounds on the NSI parameters |e„ r | and 


6 T t £ 


A iT I 


in the framework of the so-called two-flavor hybrid mode 


SK I and SK II data sets, the 90 % C.L. upper bounds are given by 


3 , 


rom the combined 


251 


C /XT | 


< 0.033, 


I err - < 0.147. 


( 4 ) 


The MINOS experiment has also used its data to set the following bound —0.2 < e MT < 0.07 


at 90 % C.L. 


27]. However, this bound is weaker than the one set by SK. 


In what follows, we will work with the exact three-flavor neutrino oscillation probabilities 
and consider the following ranges for the relevant NSI parameters: 3 


-0.3 < e eM < 0.3, 
—0.5 < e er < 0.5 , 
-0.04 < e^r < 0.04, 
-0.15 < e TT < 0.15, 


( 5 ) 


3 We choose a smaller range than the current 90 % C.L. allowed range in Eq. ([3]), since outside these ranges 
the x 2 corresponding to the ICAL data is already large. 


5 











FIG. 1: The relative probability difference as a function of the NSI parameters e^ T and 

e TT for cos 6 = —1 (left panel) and cos 9 = —0.55 (right panel). Please note the scale of the 
colorbars to the right of each panel. The following values of the neutrino parameters have been 
used: 9 12 = 34°, 613 = 9.2°, #23 = 45°, <5 = 0 (no leptonic CP-violation), Ara^ = 7.5 • 10~ 5 eV 2 , 
and Ami = +2.4 • 10 ~ 3 eV 2 (normal neutrino mass hierachy). All NSI parameters, except e, tT and 
e TT , have been set to zero. 


while for the oscillation parameters we assume the following true values: 

sin 2 0 1 2 = 0.31, sin 2 6*23 = 0.5 , sin 2 20 13 = 0.1, 

Amh = 7.5 x 1CT 5 eV 2 , lAmJjJ = 2.4 x 1CT 3 eV 2 , <5cp = 0 . ( 6 ) 

The NSI parameter(s), if present, will alter the neutrino oscillation probabilities. Oscil¬ 
lograms showing the impact of NSI on the neutrino oscillation probabilities have appeared 
vastly in the literature. The muon neutrino survival probability P w , is affected most by the 
parameters | e /tr | and |e TT — e MAt |, while the transition probability P efl depends on |e e/i | and 
|e Mr |. This dependence can be used to discover NSI parameters using neutrino oscillation 
data or constrain them. 

It is also known that the dependence of the neutrino oscillation probabilities on NSI 
parameters is different for the NH and IH cases. Indeed, since measurement of the neutrino 
mass hierarchy is one of the prime goals of the INO experiment, it is pertinent to ask 
how the mass hierarchy sensitivity of the experiment alters in the presence of NSI. For the 
mass hierarchy determination, what matters is the difference in the oscillation probabilities 
between NH and IH. Therefore, it is appropriate to ask how the difference in the probabilities 
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FIG. 2: The relative probability difference as a function of the NSI parameters e efl and e er for 
cos 9 = — 1 (left panel) and cos 6 = —0.55 (right panel). The neutrino parameter values used are 
the same as in Fig. |T] except that e efl and e eT are non-zero, while all other NSI parameter values 
have been set to zero. 


between NH and IH changes in presence of NSI. In order to show the impact of NSI on the 
mass hierarchy sensitivity, we present in Figs. [T] arid [2] the contours of the relative difference 
A^ between the neutrino oscillation probabilities P Q( g (including NSI) corresponding to NH 
and IH. We define the relative probability difference A 7 as follows (c/. the definition of the 
total CP-asymmetry in Ref. 28]) 

A-P«/(fy _ /AF™(£.9JdB 


dMH 

■tt-af3 


(*.) = 


/£p a 7(p,0 2 )dp’ 


(7) 


where 


ap7(p, e z ) = pJ(E, e z ) - p™(E, e z ), 
sp 7 ( p , 0 ,) = p^{E,e z ) + PS(E,e z ), 


P^ and P™ being the u a —> up oscillation probability for NH and IH, respectively. In 
each case, we calculate A 7 for a particular zenith angle 9 Z , while the energy dependence 
is integrated out in the range (1-100) GeV. 

In Fig. m we show the relative probability difference H])]] 1 in the e MT -e TT plane, keeping 
e e/i and e er fixed at zero. 4 The probabilities are calculated numerically within the full three- 


4 Here and throughout the rest of this work, we keep = 0. 
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generation oscillation paradigm, using the PREM 29] density profile for the Earth matter. 
We compute this for two benchmark zenith angles of cos 0 2 = —1 and —0.55 corresponding 
to neutrino baseline lengths of L — 12742 km and L = 7000 km, respectively. The colors 
represent the contours corresponding to the values of shown in the colorbar. The (0,0) 
point in the e MT -e TT plane corresponds to neutrino oscillations without NSI (i.e. standard 
neutrino oscillations). At all other points, NSI are included in the model, and this can be 
observed to alter the mass hierarchy sensitivity parameter A^ H . Note that for standard 
oscillations we have A^ ~ —5%. This small relative difference is what the atmospheric 
neutrino experiments observing are exploiting to determine the neutrino mass hierarchy. 
When the NSI parameters are switched on, A^ H changes. The relative difference A^ is 
seen to increase for e MT < 0 and decreases further to larger negative values for e MT > 0 . 
However, since for the hierarchy measurement what is relevant is the absolute difference 
I^TI , we will later see that for all | e /tr | > 0 , the hierarchy sensitivity increases as long as all 


the parameters are kept fixed between the NH and IH cases. The change in Jy with |e 




is seen to be in the same direction for both the zenith angles shown, though its magnitude 
is seen to be larger for cos 9 Z = — 1 case. The change in with e TT is less pronounced. 
In particular, for the cos 9 Z = —0.55 case, the dependence on e TT is marginal. Even for the 
cos 6*2 = —1 case, the dependence of on e TT for e TT = 0 is negligible. For larger values 
of |e Tr |, the role of |e TT | is to reduce the overall change in due to NSI, and this happens 
for both positive and negative e rr . 

In order to understand the change of the probability difference as a function of the NSI 
parameters, we can series expand the neutrino oscillation probabilities in orders of the NSI 
parameters and keep only the first-order terms. The expression for the difference in the 
muon neutrino survival probability between NH and IH, keeping only leadin g-or der terms 


in NSI parameters and neglecting the standard matter effects, is given by 

A 


30 


AP^ H ~ -2Re(e MT ) sin 20 23 f sin 2 20 23 ^ sin 


| I L 


6 31 1 

2 E 


+ 4 COS 2023 


| Am 


sm 


32] 

2 I Am: 


2 I 
311 


4 E 


+ ( | W | - I'£rr |) Sin- 2023 COS 2023 


f AL . | Amp 1 |L 


sm 


-4- 


A 


sm - 


I Am!,| L 


■ ( 8 ) 


\2E 2 E lAmpJ 4 E 

where A = 2y/2GF?i e E. Note that Eq. (EJ) depends only on the parameters Re(e^ T ) and 
|e wt | — |e TT | to leading order. Dependence on e efl and e eT appear only at higher orders, which 
can be neglected unless these parameters are taken to be large. Therefore, in Fig. [Q we 
show the relative probability difference for the survival channel in the e /tr -e rr plane keeping 












the other NSI parameters at zero. For e TT = 0, the expression clearly shows that |AP^ H |, 
and hence |A^|, grows with |e MT | and flips sign when the sign of e tlT changes. The quantity 
A^ H is positive for e Mr < 0 and negative for e /tr > 0. This agrees fairly well with the exact 
results shown in Fig. |T] The impact of e TT on the other hand is more involved. Using Eq. (JSJ), 
we note that for any given large value of e )lT , we should have the highest possible |A^ H | for 
e TT = 0, and since the dependence on this parameter comes in the form of |e Tr |, we should 
have lower |A^| on both sides of e TT = 0. On the other hand, for e^ T = 0, the AP^ H 
obtains contribution only from the second term, and there is a relative sign between the two 
terms in the parentheses. As a result for e /tT = 0 we do not expect large contribution to 
AP“ H from NSI. These features can be observed in the exact result in Fig. [Tl 

In Fig. [21 we present the A^ contours in the e e ^-e eT plane with e MT and e TT fixed at 
zero. The probability P eM depends crucially on the NSI parameters e efl and e eT at leading 
order, and hence, NSI bring significant change to |A^ H |. In this case, the corresponding 
analytic expression is complicated and we refer the reader to Ref. 0 for a related expres¬ 
sion for the approximate formula. However, the exact results shown in the figure tell us 
that the presence of the NSI parameters e e/i and e eT could bring substantial change to the 
mass hierarchy sensitivity of atmospheric neutrino experiments. In fact, e e/1 and e eT could 
either increase or decrease the mass hierarchy sensitivity compared to what we expect from 
standard oscillations. 

III. EVENT SPECTRUM AT INO 

The ICAL (Iron CALorimeter) detector at INO will be a 50 kton detector with layers 
of magnetised iron interleaved with glass Resistive Plate Chambers (RPC), which will serve 
as the active detector element. The atmospheric neutrinos in z/ M , zyf, u e , and W e species 
come from decay of pions and kaons produced from cosmic ray interactions with the Earth’s 
atmosphere. These neutrinos can interact with the detector nucleons producing the corre¬ 
sponding charged lepton. The dense iron material of ICAL helps to detect muons through 
their long tracks 5 , while the magnetic held allows the identification of their charge. Since 
Earth matter effects develop only in either the neutrino or the antineutrino channel for a 

5 The electrons give rise to an electromagnetic shower in the detector, which cannot travel far and is 
therefore more difficult to observe in this class of detectors. 
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given mass hierarchy, this charge identification capability gives ICAL an edge to better ob¬ 
serve the Earth matter effects, and hence, the neutrino mass hierarchy. The capability of 
this experiment to help discover the mass hierarchy has been studied before by the INO 
collaboration using information on muon energy and angle in Ref. Q and using both the 
muon energy and angle information as well as hadron energy information in Ref. js] . In 
this work, we only use the muon energy and angle information and quantify the change in 
the mass hierarchy sensitivity of ICAL in presence of NSI. We also study the prospects of 
constraining or discovering the NSI parameters with the muon event sample. 

For calculating the predicted number of /i - and /i + events in ICAL, we use the same pre¬ 
scription as in the earlier INO collaboration papers. The unoscillated events are calculated 
using the Nuance event generator modified for ICAL. The oscillation probabilities, with and 
without NSI effects, are implemented through a re-weighting algorithm. Finally, the muon 
reconstruction efficiency, charge identification efficiency, and muon energy and angular res¬ 
olutions are folded in as described in Refs. [ 2 -6]. The new ingredient in the simulations 
performed for this work is that while all the earlier papers used muon sample in the energy 
range (1-11) GeV, we consider in this work a much larger energy range of (1-100) GeV. In 


order to do that, we extend the earlier study 


33] for muon detector response to 100 GeV 


from detector simulations done with the GEANT-based code developed for ICAL. The muon 
energy and zenith angle resolutions, as well as the charge identification efficiency and recon¬ 
struction efficiency are obtained as a function of muon energy and zenith angle. This is then 
folded with the oscillated events to obtain the final muon spectrum expected in ICAL. We 
generate raw events corresponding to 1000 years of running of ICAL in order to reduce the 
Monte Carlo fluctuations and normalise the final events to 10 years of running. This event 
sample is then binned in energy and zenith angle bins as follows. For the zenith angle we 
have 20 equal size bins in the cos6L range (—1,1). For the energy, we take variable bin sizes 
to ensure that there are reasonable number of events in each bin. Between muon energies 
(1-10) GeV, we take 9 energy bins of size 1 GeV, and between (10-100) GeV, we take 3 
energy bins of size 30 GeV. 

The predicted events are then analysed by a statistical procedure identical to the one 
used in the earlier papers. A y 2 function is defined as 


X=X 2 (/0 + X> + ), 


(9) 
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Ni Nj r 


x 2 (^) = 


c± ■“—' ■“—' 

i=\ j =1 






/ U ^ 

V JVg-CA> 


E^(io) 


fc=l 


«?>*) = iv'*‘(/< ± )(i + E>«£) + °K. 

' fe=l ' 


± 2 ] 
fc J 


( 11 ) 


N'ij (l-i) and A^|j x (/i ± ) being the predicted and ‘observed’ number of p, events in ICAL, 
respectively, the correction factors due to the k th systematic uncertainty, and the 
corresponding pull parameters. In this analysis, we include 5 systematic uncertainties. These 
are 20 % error on flux normalisation, 10 % error on cross-section, 5 % uncorrelated error 
on the zenith angle distribution of atmospheric neutrino fluxes, 5 % tilt error, and a 5 % 
overall error to account for detector systematics 6 . The individual contributions from fi~ and 
/i + data samples are calculated by minimising Eq. (TTOT) over the pull parameters. These are 
then added [cf. Eq. (IIUI) ] to obtain the y 2 for a given set of oscillation (and NSI) parameters. 
This resultant y 2 is then marginalised over the oscillation parameters, and when applicable, 
over the NSI parameters. We assume for the oscillation parameters the true values given in 
Eq. d6j) and marginalise our y 2 over their current 3cr ranges. We include priors defined as 

x ^ = ( ? ^r 1 ) 2 ’ (12) 

where p t me is the assumed true value of the parameter p and a p the la error on it. We include 
oriors on sin 2 0 2 3 ; and sin 2 20 13 with la errors of 1 %, 2 %, and 0.005, respectively 


M 


35]. For the NSI parameters, the y 2 is marginalised over their range given in Eq. 


For all results presented in this work, we use 500 kton-year of statistics in ICAL. We next 
define the different y 2 that we compute in this work for the different physics studies we 
perform. 

Sensitivity to Neutrino Mass Hierarchy: To find the sensitivity of ICAL to the 

neutrino mass hierarchy, we compute the y 2 according to Eqs. ([9]) and ffTOlh where 
correspond to the right hierarchy and A r '^( / u ± ) correspond to the wrong hierarchy. The 
Ay^n thus obtained is then marginalised over the oscillation and NSI parameters. We do 
this for different assumed true values of the NSI parameters. 


6 Simulations to estimate the detector systematic uncertainties in ICAL is underway. This number could 
therefore change when better estimates of this become available. 
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Bounds on NSI parameters: In the event that there is no signal for NSI in the ICAL 
data, one will be able to give upper bounds on the NSI parameters at a given C.L. In order 
to find the expected sensitivity of ICAL to constrain the NSI parameters, we compute the 
y 2 by generating for standard oscillations (with all NSI parameters set to zero) and 

fitting this with (/i ± ) computed with non-zero NSI parameters. The corresponding Ayg, 
where S stands for sensitivity, obtained after marginalising over oscillation gives a measure 
of the sensitivity reach of ICAL to NSI. 

Discovery of NSI parameters: If on the other hand, one finds a signal of NSI in 
the ICAL data, this would be a discovery of NSI, and hence, physics beyond the Standard 
Model. Of course, the NSI parameters in this case have to be above a certain value to be able 
to produce a discoverable signal at ICAL. We find this lower limit on the NSI parameters 
needed to be discovered at ICAL for a given C.L. by generating with NSI and 

fitting this with corresponding to standard oscillations. The corresponding Ay^ 

obtained after marginalising over oscillation gives a measure of the discovery reach of ICAL 
to NSI. 

Precision on NSI parameters: Finally, for a given set of NSI parameters, one can use 
the ICAL data to produce C.L. contours in the NSI parameter space. We will show these 
contours at the 68 %, 95 %, and 99 % C.L. in the e MT -e TT and e e ^-e eT planes. For this, we 
will generate A®* (/.A) for a certain set of NSI parameters and fit it with all values of the 
NSI parameters in the given plane while marginalising over the oscillation parameters. 

IV. EXPECTED RESULTS FROM INO 

We now present our main results. We first show the impact of NSI on the mass hierarchy 
sensitivity of ICAL, which is the main goal of the experiment. We next give the sensitivity 
reach of this experiment in constraining NSI parameters. Subsequently, we quantify the 
NSI discovery potential at INO. Finally, we briefly discuss with what precision the NSI 
parameters could be measured at INO if they were indeed above the discovery limit. 
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FIG. 3: The Ay^j H , giving the expected mass hierarchy sensitivity from 10 years of running of 
ICAL, as a function of the true value of NSI parameters. We keep only one e Q( a(true) to be non-zero 
at a time, while others are set to zero. The Ay 2 is obtained as explained in the text. However, the 
resultant Ay 2 is not marginalised over the oscillation parameters as well as NSI parameters. 


A. Impact of NSI on Mass Hierarchy Sensitivity 


We noted in Figs. [Hand[2] that the difference in the oscillation probabilities between NH 
and IH changes in the presence of NSI. Since the mass hierarchy sensitivity is defined in 
terms of the difference in the event distribution for NH and IH, it is therefore expected that 
the mass hierarchy sensitivity of the experiment would change in the presence of NSI. The 
mass hierarchy sensitivity for standard oscillations using only the muon events in ICAL is 
given in Ref. [2]. We revisit the mass hierarchy sensitivity in ICAL in the presence of NSI 
parameters and show our results in Figs. [3] and [4[ The Ay 2 m corresponds to the difference 
in y 2 of the fit with the wrong and the right hierarchy as a function of the true value of the 
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FIG. 4: The Ay^j H , giving the expected mass hierarchy sensitivity from 10 years of running of 
ICAL, as a function of the true value of NSI parameters. We keep only one e aj a(true) to be non¬ 
zero at a time, while others are set to zero. The Ay 2 is obtained after marginalisation over the 
oscillation parameters as well as NSI parameters as explained in the text. 


NSI parameter. For the sake of simplicity, we take only one non-zero NSI parameter in the 
data at a time. For instance, the black curves in the top-left panels of Figs. [3] and [4] are 
obtained as follows. The data are generated for NH and a given true value of e eM (shown 
as the x-axis). The oscillation parameters in data are taken from Eq. (E]) and all other NSI 
parameters are set to zero. This is then fitted with a theory prediction corresponding to IH. 
In Fig. [3l we present the Ay^ H obtained when all oscillation and NSI parameters in the fit 
are fixed at their assumed true values. In Fig. [Tj we marginalise the Ay 2 over the oscillation 
parameters lAmjjJ, sin 2 023, and sin 2 2013 with priors. The Ay 2 IH is also marginalised over 
the NSI parameter which is taken as non-zero in the data, while the other NSI parameters 
are kept fixed at zero. For instance, in the top-left panel, the Ay 2 is marginalised over e efl , 
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while the other NSI parameters are kept fixed at zero. In all cases, the Ay 2 IH is marginalised 
over the oscillation parameters lAm^l, sin 2 9 2 3 , and sin 2 0 13 with priors included as described 
in the previous section. The other panels are also obtained in a similar way. 

The horizontal black dashed lines in the four panels of Figs. [3] and |4] show the mass 
hierarchy sensitivity expected in ICAL for the case when there are no NSI considered in 
either the data or the fit. A comparison of this with the black solid and red dashed curves in 
the figure reveals that presence of NSI in the data could change the mass hierarchy sensitivity 
of ICAL. I 11 particular, we see that the Ay^ H changes sharply with the true value of e efl and 
e eT . In presence of NSI, we note that the Axl/m increases for e e/t (true) > 0 and e eT (true) > 0, 
while it decreases for e eAl (true) < 0 and e eT (true) < 0 , compared to what is expected for 
standard oscillations. 

These features can be understood from Figs. |Tjand[2] Since the NSI parameters e eM and 
e er mainly affect the appearance channel P efl , we refer to Fig. [2] to understand the upper 
panels of Figs. [3] and [4] The (0,0) point of Fig. prefers to standard oscillations and gives 
the mass hierarchy sensitivity shown by the black dashed lines in Figs. [3] and [4l If we stay 
on e er = 0 and change e e/1 , we note from the left panel (cos 6 Z = —0.55) of Fig. [2] that for 
e efl 0 |A^ H | decreases, while for e eil > 0 it increases. This is less clear in the core-crossing 
bin, however, since the largest mass hierarchy sensitivity at ICAL comes from zenith angle 
bins close to cos 9 Z = —0.55, this feature stays in the final Ax^ H . 

The effect of e TT on the mass hierarchy sensitivity is seen to be less severe from the bottom- 
right panels of Figs. [3] and [U This NSI parameter affects the muon neutrino survival channel 
the most. Figure [Q reveals that the impact of e TT on |A^| (when e MT = 0) is very small for 
both the core-crossing and the cos 9 Z = —0.55 bin. The impact of e lT on the mass hierarchy 


sensitivity is more interesting and has been discussed in Refs. 10-12], For e TT = 0, we see 
that |A^ H | could change up to 20 % for the core-crossing bin and a few percent for the 
cos 6*2 = —0.55 bin, as we change e MT . Note that for standard oscillation |A^| is already a 
very small number, and hence, the relative change of |A^| due to e flT is significant. This is 
reflected in the bottom-left panel of Fig. [3] where we see a large increase in Ay^[ H with e MT . 
However, once we marginalise over oscillation and e MT in the fit, this increase is washed out 
and we obtain no significant impact of e MT on Ay^ H in the bottom-left panel of Fig. [H 
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FIG. 5: The A%|, giving the sensitivity reach of 10 years of ICAL data in constraining the NSI 
parameters in the event that the data show no signal of any new physics, as a function of the fit 
value of NSI parameters. We keep only one e aj g(true) to be non-zero at a time, while others are 
set to zero. The has been marginalised over the oscillation parameters and NSI parameters 
as explained in the text. 


B. Expected Bounds on NSI 

In Fig. El we present the expected sensitivity reach of ICAL in constraining NSI param¬ 
eters. The figure shows the expected sensitivity for each of the NSI parameters e efl (top-left 
panel), e er (top-right panel), e Mr (bottom-left panel), and e TT (bottom-right panel). This 
figure is obtained as follows. We use as data the event distribution at ICAL corresponding 
to standard oscillations by setting all NSI parameters to zero. This is then fitted with the 
predicted event distribution which includes one non-zero NSI parameter at a time, and the 
corresponding Ay§ calculated. In Fig El we show this A\g as a function of the NSI param¬ 
eter that is allowed to be non-zero in the fit. The Ay| is marginalised over the oscillation 
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FIG. 6: The improvement in the expected bounds on NSI from increasing the considered muon 
energy range from 20 GeV to 100 GeV in the analysis. We show this only for the NSI parameters 

€[1T ■ 


parameters sin 2 # 23 , and sin 2 26 Priors on the three oscillation parameters were 

included as described in the previous section. The resultant Ax§ shows little change as a 
result of marginalisation over them. The black solid curves are obtained when the data are 
considered corresponding to NH, while the red dashed curves are for data corresponding to 
IH. We keep the hierarchy fixed to its assumed true value in the fit. 

The expected sensitivity for IH is only marginally worse than that for NH. At the 90 % 
(3cr) C.L., the expected bounds on the NSI parameters from 500 kton-years of statistics in 
ICAL for NH can be read from the figure as 


- 0.119 (-0.3) < e c/1 < 0.102 (0.2), 
-0.127 (-0.27) < e eT < 0.1 (0.23), 
-0.015 (-0.027) < e^ T < 0.015 (0.027), 
-0.073 (-0.109) < e TT < 0.073 (0.109). 


For the IH case, the bounds are comparable and can be read from the figure. 

In Fig. [6] we show the improvement that we obtain in the sensitivity reach of ICAL to 
the NSI parameters when we increase the muon energy range considered in the analysis from 
20 GeV (red dashed curve) to 100 GeV (black solid curve). The 3a bound on e MT improves 
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FIG. 7: The Ay^, giving the discovery potential of 10 years of ICAL data in observing NSI, as a 
function of the true value of NSI parameters. We keep only one NSI parameter to be non-zero at a 
time, while others are set to zero. The Ay^ has been marginalised over the oscillation parameters 
as explained in the text. 


from —0.033 < e MT < 0.033 to —0.027 < e MT < 0.027, when we increase the muon energy 
from 20 GeV to 100 GeV in the data. 

C. Discovery Reach for NSI Parameters 

In the previous section, we looked at how well ICAL will be able to constrain NSI pa¬ 
rameters if its data were consistent with just standard oscillations. Next, we take the com¬ 
plementary view and ask ourselves that if NSI parameters were indeed non-zero, at what 
C.L. would ICAL be able to tell them apart from standard oscillations. In other words, 
we are looking for the limiting true values of the NSI parameters above which the data at 
ICAL would be a signal for NSI at a certain C.L. For that, we now consider data for various 


18 












































(assumed) true values of the NSI parameters and fit them with a predicted event spectrum 
corresponding to standard oscillations and compute the resultant Ayjy We present this in 
Fig. [3 For simplicity, we consider only one non-zero NSI parameter at a time in the data. 
We marginalise over the oscillation parameters IA 777 . 3 -J, sin 2 0 2 3 ; and sin 2 # 13 with priors im¬ 
posed on each one of them as discussed before. The black solid curves correspond to the 
case for NH, while the red dashed curves are for IH. We keep the hierarchy to be the same 
in the theory as in the data. The figure shows the expected Ay^ for the discovery of each of 
the NSI parameters e e/i (top-left panel), e eT (top-right panel), e ^ (bottom-left panel), and 
e TT (bottom-right panel). While the nature of the curves are very similar to the ones we 
had in Fig. 0 the values of the Ay^ are different. With 500 kton-years of data, the ICAL 
experiment will be able to give a signal of NSI at the 90 % (3cr) C.L. for NH if 


e e#1 < -0.116 (-0.28), 
e er < -0.12 (-0.29), 
e MT < -0.015 (-0.027), 
e TT < -0.07 (-0.104), 


e e/i > 0.105 (0.2), 
e eT > 0.102 (0.23), 
e Mr > 0.015 (0.028), 
e TT > 0.07 (0.103). 


The corresponding limiting values for IH are similar, as can be seen from the figure. 


D. Precision on NSI Parameters 

In Figs. [8] and [9l we show the projected C.L. area in the NSI parameter space allowed 
after 10 years of running of the ICAL experiment. In Fig. El we show these for the case where 
we assume that there are no NSI, or in other words, when the assumed true values of the 
NSI parameters are taken as zero, shown by the black dots in the figure. The black dotted, 
blue dashed, and red solid contours show the 68 %, 95 % and 99 % C.L. in two-dimensional 
NSI planes. The contours are marginalised over the oscillation parameters after including 
priors that are described above. The NSI parameters other than the ones appearing in the 
two-dimensional plane are kept fixed at zero. NH is assumed for all plots. The corresponding 
contours for the IH case are very similar and we do not repeat them for brevity. 

In Fig. El we show the C.L. contours for the case where we assume the true values of 
the NSI parameters to be non-zero. The upper panels show the C.L. in the e^-err plane 
when the true values of (6^,6^) are taken as (0.02,0), (0,0.075), and (0.02,0.075) for the 
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FIG. 8: The expected C.L. contours in the given NSI parameter plane. The other NSI parameters 


are set to zero. The NH has been assumed to be true. The black dots show the points where the 


data were generated, which are for no NSI in this case. 










FIG. 9: The expected C.L. contours in the given NSI parameter plane. The other NSI parameters 
are set to zero. The upper panels are drawn in the e^ r — e TT plane, while the lower panels are 
drawn in the e er — e e/Jj plane. The black dots show the points where the data were generated. 
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left, middle, and right panels, respectively. These assumed true points are shown by black 
dots in the plots. The values e efl and e eT are assumed to be zero in both the data as well as 
in the fit, and we do not show the contours in planes involving these parameters. The lower 
panels are similar to the upper panels except that now we show the C.L. in the e eT -e eil plane 
when the true values of ( e eT ,e efl ) are taken as (0,0.1), (0.15,0), and (0.15,0.1) for the left, 
middle, and right panels, respectively. For these panels, the values e MT and e TT are assumed 
to be zero in both the data as well as in the fit. Again, the figures are for NH, however, the 
ones for 1H are similar and we do not present them for brevity. 


V. CONCLUSIONS 


The study of the physics potential of the ICAL detector at the proposed India-based 
Neutrino Observatory is underway. As a part of this on-going effort, we probe in this work 
the impact of NSI parameters on the expected signal at ICAL and the physics conclusions 
that one can draw out of it. The neutral-current NSI if present, could alter the propagation 
of atmospheric neutrinos inside the Earth matter changing the signal at ICAL. This change 
due to NSI can be used to study the NSI parameters. On the other hand, one needs to 
estimate how much the potential of ICAL to standard physics gets modified in the presence 
of NSI. In this work, we have taken both these considerations into account and studied the 
physics potential of ICAL in presence of NSI. 

Measurement of the neutrino mass hierarchy is the primary goal of the ICAL atmospheric 
neutrino experiment. We showed how the difference in the neutrino oscillation probabilities 
between NH and IH change in presence of NSI. We defined the relative probability difference 
A for the oscillation channel u a —> and showed the oscillograms for A ^ and A^ H , the 
two oscillation channels relevant for the atmospheric neutrinos in ICAL. These oscillograms 
(and all other results shown in this paper) were obtained from an exact numerical calculation 
of the three-generation neutrino oscillation probabilities using the PREM profile for the 
Earth matter density 29]. The oscillograms show that the relative difference A ^ changes 
significantly with e MT compared to its Standard Model value, while A^ t H is seen to vary 
sharply with the values of e eM and e eT . The impact of the NSI parameter e TT is seen to be 
less important. 

We next simulated fi~ and /i + events in the ICAL detector in presence of NSI and defined a 
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X 2 function, including energy and zenith angle correlated as well as uncorrelated systematic 
uncertainties, to give C.L. predictions for the estimated sensitivity of ICAL.' The y 2 is 
marginalised over the NSI parameters and the oscillation parameters lAmjjjJ, sin 2 # 23 , and 
sin 2 20 13 after putting priors on them. Using this we presented the change in Ay 2 IH if NSI was 
assumed to be a certain true value in Nature. We showed that the Ayjq H increases rapidly 
for e e/t (true)> 0 and e eT (true)> 0, while it decreases for e e/i (true)< 0 and e eT (true)< 0 
compared to what we expect for standard oscillations. This behavior can be understood 
from the oscillograms we showed. The impact of the NSI parameter e TT is small, however, 
the Ay^ H could vary significantly with e MT . However, if we allow for marginalisation over 
the oscillation (especially lAmg-J) and NSI parameters, the Ay^ H comes to be around the 
value predicted by the Standard Model. 

We next showed the potential of ICAL in discovering or constraining NSI. If the case that 
ICAL was consistent with no NSI in the data, we presented the expected upper limit on the 
NSI parameters. At the 90 % (3cr) C.L. we have for the NH the limits 


-0.119 (-0.3) 
-0.127 (-0.27) 
-0.015 (-0.027) 
-0.073 (-0.109) 


< e efl < 0.102 (0.2), 

< e er < 0.1 (0.23), 

< < 0.015 (0.027), 

< e rT < 0.073 (0.109). 


The limits for IH are similar. Compared to the current 90% C.L. bounds given in Eqs. ([3]) 
and (|4]) the expected bounds from ICAL are promising. We next considered the case where 
the data at ICAL is consistent with NSI and we gave the expected statistical significance 
with which ICAL will rule out the theory with no NSI. We calculated the range of the 
NSI parameters that would lead to 90% (3a) C.L. discovery of NSI at ICAL. Finally, we 
presented the C.L. contours in the two-parameter NSI planes, for different choices of true 
values of the NSI parameters. 


7 It has been shown that the inclusion of hadron energy information in the analysis of ICAL data improves 
the mass hierarchy sensitivity of ICAL. The impact of the hadron energy information on the sensitivity 
of ICAL to NSI is being studied in an independent work [20 1. 
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